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1. Introduction: The hypothetical Hilbert-Polya operator 

There is an old idea, usually attributed to Hilbert iJj and Polya [4J that the nontrivial 
(i.e. complex) zeros s„ of the Riemann zeta function ({s) have a spectral interpretation. 
Writing s„ := | — ir„, the Riemann hypothesis states that the nonimaginary solutions 
t:„ of C(i — ir„) = are real, that is the nontrivial zeros s„ lie on the critical line 
Res = |. The Hilbert-Polya approach towards a proof of the Riemann hypothesis 
consists in finding a Hilbert space Ti and a self-adjoint operator H inTi whose discrete 
spectrum is exactly given by the nontrivial zeros r„ = i (s„ — |). 

Around 1950, Selberg [5] introduced his zeta function Z{s) in analogy with ({s) and 
with the intention to shed some light on the nontrivial Riemann zeros and the Riemann 
hypothesis. He noticed the striking similarity between his famous trace formula for 
the Laplace-Beltrami operator on e.g. compact Riemannian manifolds and the explicit 
formulae of number theory, whose most general form is Weil's explicit formula [3]. 
The nontrivial zeros of the Selberg zeta function Z{s) fulfil the analogue of Riemann's 
hypothesis and appear in the spectral side of the trace formula being directly related 
to the spectrum of the Laplacian. The other side of the trace formula has a purely 
geometrical interpretation, since it is given by a sum over the length spectrum of the 
closed geodesies (periodic orbits) of the geodesic fiow, i.e. the free motion of a point 
particle on a given hyperbolic manifold. This system was already studied by Hadamard 
[6l [7| in 1898 and has played an important role in the development of ergodic theory 
ever since. Hadamard proved that all trajectories in this system are unstable and 
that neighbouring trajectories diverge in time at an exponential rate, the most striking 
property of deterministic chaos. 

In 1980, Gutzwiller |8] drew attention to this system as a prototype example 
of quantum chaos by identifying the Laplacian on hyperbolic manifolds with the 
Schrodinger operator in quantum mechanics. In this way he related the nontrivial 
zeros of the Selberg zeta function to the quantum energies of a dynamical system whose 
classical trajectories are chaotic. Furthermore, he realized that the Selberg trace formula 
is an exact version of his trace formula, the celebrated Gutzwiller trace formula [9j, which 
holds for general quantum systems with a chaotic classical counterpart, but in this case 
only approximately, i.e. in the so-called semiclassical limit where Planck's constant h 
approaches zero. 

In 1985, Berry [10] emphasized that the search for the hypothetical Hilbert-Polya 
operator in terms of a Schrodinger operator obtained from the quantization of a 
classically chaotic system might be a fruitful route to proving the Riemann hypothesis. 
He discussed in detail the properties of this operator that are suggested by the quantum 
analogy. Prompted by a paper written by Connes [llj (see also [12]), who devised 
a self-adjoint operator (Perron-Frobenius) of a classical dynamical system together 
with a classical (Lefschetz) trace formula in noncommutative geometry. Berry and 
Keating [T3l [T4] speculated that the conjectured Hilbert-Polya operator might be some 
quantization of the extraordinarily simple classical Hamiltonian function Hci{x,p) of a 
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single coordinate x and its conjugate momentum p: 

Hdix,p) := xp. (1) 

Inspired by [T5lfT6] . Berry and Keating [1^ suggested to investigate quantum graph 
models of ([1]), in particular the spectrum of these operators. One of the first researchers 
who dealt with differential operators on graphs was Roth [T7] who derived a trace formula 
for the heat kernel of the Laplacian with Kirchhoff boundary conditions. Von Below 
|18] considered the heat equation on graphs and derived a characteristic equation for the 
eigenvalues of the weighted Laplacian on graphs. Some physical quantum graph models 
were considered by Exner and Seba [19] who discussed i.a. the scattering problem for a 
free quantum particle on a star graph. A method to approximate mesoscopic systems 
like thin branching systems by quantum graphs was discussed by Exner and Post [20] 
and Post |2T]. Carlson |22] used semigroups on graphs to simulate the blood flow in the 
human arterial system. Kottos ans Smilansky [T5| [T6] introduced quantum graphs as a 
model for quantum chaos. 

In this paper, we study the quantization of the classical Berry-Keating Hamiltonian 
([1]) in the Hilbert space L^(]R>, dx) and on compact quantum graphs and give a complete 
classification of the self-adjoint realizations of the corresponding Berry-Keating operator. 
In addition, we also study the quantization of the corresponding "squared" operator. It 
turns out that no self-adjoint realization of ([1]) exists which yields as eigenvalues the 
Riemann zeros. 



2. Classical dynamics and quantization of the Berry-Keating operator 

Let us consider the classical dynamics of a particle moving on the real line R generated by 
the Berry-Keating Hamiltonian ([1]) with corresponding phase space P : (x,j9) G M x R. 
The classical time evolution (Hamiltonian flow) is governed by Hamilton's equations 

i(t) = ^ = x(t) and p[t) = -^ = -pit). (2) 

Starting at time t = at an arbitrary point {xo,Pq) G P in phase space, the unique 
solutions are [13j 

x(t) = Xqc* and p(t) = Poe~^. (3) 

Obviously, the point (0, 0) G P is an unstable point. We note that the Hamiltonian 
([1]) is time independent corresponding to the conserved "energy" E := Hci{x{t),p{t)) = 
XqPo G R, and thus the particle moves in V on the "energy surface" (hyperbola) xp = E. 
Obviously, the classical motion is unbounded. Therefore, Berry and Keating [T3| [H] 
introduced some regularization procedures, leading to a truncation of phase space, which 
we shall discuss below, but first we would like to discuss quantum mechanics. 

Quantization of the classical system requires to choose a Hilbert space "H and to 
replace the classical Hamiltonian (II]) by a self-adjoint operator H in l-i. With the 
standard choice l-L := L^(R, dx), the simplest operator corresponding to ([1]) is obtained 
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by Weyl ordering of the coordinate operator x (acting by multiphcation) and the 

''da; 



momentum operator p = — i^4- (acting by differentiation) leading to the Berry-Keating 



operator [131 E] 



^BK ■=^{xp + px) = -ih + ^ , (4) 

and the Schrodinger equation 

ih^^^^ = H^K^ix,t). (5) 

As was to be expected from our discussion of the classical motion, the operator -f^BK is 
unbounded and does not have a discrete spectrum corresponding to bound states, but 
rather has a continuous spectrum A G R corresponding to scattering states obtained by 
solving the eigenvalue problem 

HsK'ipix) = Xipix). (6) 

Writing A = hk, A; G M, Planck's constant drops out from ([6]), and the eigenvalue 
problem reads (s := — | + ifc) 

= s^.(,). (7) 

For a; G M, (E]) possesses the general solution 

0,(x) = CiX+ + C2Xi, (8) 
where denote the generalized functions (see e.g. |23| p. 87]) 

{0 for X < f Ixl"* for X < 

and xl := < , (9) 

x'* for X > [0 for X > 

which is well defined for Re s > —1. In [13j, Berry and Keating studied as a special case 
the simplest choice for the continuation of the eigenfunctions across the singularity at 
X = by considering the even eigenfunctions (ci = C2 = c) 0g™'^(x) = c|x|*. 

Let us discuss in more detail the case that the quantum dynamics takes place on the 
positive half-line x G ]R>. Then Hbk acting on D(R>), the set of infinitely continuous 
differentiable functions with compact support on M>, is essentially self-adjoint (see e.g. 
[24j [both deficiency indices are equal to zero]). Therefore, the closure of this operator 
is self-adjoint. The general solution of the time-independent Schrodinger equation ([6]) 
is then given by 

tpkix) := -^xl^^''" with keR, (10) 
V 27r 



which is obviously not in L^(M>,dx) and satisfies the orthonormality relation (in a 
distributional sense) 

oo 

{^k I ^fc') := J i'kix)ipk'{x)dx = 8{k - k') (11) 
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and the completeness relation 

4'k{x)'ilJi.{x')dk = 6{x — x'). (12) 
Thus, we have for any (p G L^(M>, dx) the spectral decomposition 

oo 

= j A{k)iJkix)dk (13) 

— oo 

with 

oo 

A{k) := = / ^fc(a;)0(x)da; (14) 



and (assuming (0 | 0) = 1) 

\A{k)fdk = l. (15) 



Forming a general wave packet with a given amplitude A{k) satisfying (fT3|) and ( IT5|) . 
one obtains (x G M>) 

m = y^i(lna;), (16) 

where A denotes the Fourier transform of A 

oo 

Aiy) ■■=^ J A{k)e''ydk. (17) 

— oo 

Let us also mention an alternative to the spectral decomposition f|T3|) for the Berry- 
Keating operator. Defining the Mellin transform of by 

oo 

x'-^(f){x)dx, (18) 



we obtain for the wave number amplitude A{k) (see ([Hj)) 

A{k) = ^^(l-ik], (19) 



V2 

from which 0(x) can be recovered by the inverse Mellin transform (convergent at least 
in mean square, see e.g. |25| p. 94]) 



2 

— [ 4>(s)x-'ds = — [ ^(--ik)x''^-^'''dk 
Tii J 2% J \2 J 

(20) 



5-100 
00 



A{k)iJkix)dk 
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in agreement with (fT5|) . 
The unitary group 



Uit) := expf-^i/Bx) = e-ie-*^ (21) 



h 

generated by the Berry-Keating operator (jlj) acts on functions G as (see i.a |26[ 
p. 365]) 

(f/(t)0)(x) =e-^(e~*x) . (22) 

Here we have used the relation H^k = (-^ + I)' where D := is the generator 
of scahng transformations (dilations). Let us mention that the operator D has been 
discussed by Arendt [271 [211 |29], where Ap := —D is considered as the generator of 
a semigroup on e.g. L^(M>,dx) (1 < p < oo) with Dirichlet and Neumann boundary 
conditions. 

On the other hand, the action of the unitary operator U (t) on eigenfunctions ijj of 
-f^BK gives according to 

{Um){x)=e-'i'^{x), (23) 

which in turn leads with (1221) . A = hk, s = — | + ifc and k := e^* > (t < oo) to 

iIj(kx) = kXx). (24) 

This shows that an eigenfunction ip of i^BK must be a homogeneous function with 
(complex) degree s = — | + ik. Differentiation of fl24p with respect to k, and then setting 
K = 1 leads back to the eigenvalue problem ([7]) which possesses for x G M> the unique 
solution 

For the (retarded) integral kernel KB]^{x,XQ;t) of the time-evolution operator U{t) 
one obtains {x,Xo G M>; 0(t) is the Heaviside step function) 

Kbk{x, xq] t) = e~i6 [xq - xe~*) 6(t). (25) 

We observe that the quantum mechanical time evolution follows in the configuration 
space exactly the classical trajectory (|3]). Starting at time t = with the initial wave 
function (j) G L^(M>, dx), one obtains with (123]) the wave function 4j{x, t) at a later time 
t > 



ilj{x,t) = / KBK{x,Xo;t)(j){xo)dxo 

I (26) 

= e~20 (e~*x) 

in complete agreement with (1221) . We also give the result for the resolvent kernel 
(outgoing Green's function [a small positive imaginary part (e > 0) has been added 
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to \ = hk]), see e.g. [301 p.26], 

oo 

Gbk{x,Xo;X) -.= ^1 e^(^+'^)*i^BK(x,a;o;t)dt 





oo 

^fc'(x)^fc,(xo)^^, (27) 



hk' ~X- ie 

27ri — 
n 

which satisfies the inhomogeneous time-independent Schrodinger equation (see ([6])) 

{Hbk,x - ^k) Gbk{x, Xo; hk) = 6{x -Xq). (28) 

Since the operator (jlj) acting in the Hilbert space L^(]R, dx) respectively L^(]R>, dx) 
has only a continuous spectrum, it cannot be considered (with the above realization) as 
a candidate for the hypothetical Hilbert-Polya operator. Thus, there remains the task 
to find another Hilbert space for which the quantization of the classical Hamiltonian ([1]) 
possesses a discrete spectrum. Perhaps the required space is a quantum graph, with xp 
acting on edges between vertices, a possibility already mentioned by Berry and Keating 
[13{ . It is the purpose of our paper to discuss the self-adjoint realizations on compact 
quantum graphs and in a forthcoming paper [31] we study a noncompact quantum 
graph. 

3. Semiclassical regularization of the Berry-Keating operator 

Before we come to an investigation of quantum graphs, we would like to discuss an 
alternative and very interesting approach also put forward by Berry and Keating [13] 
(see also Connes [TT| [12]) which is based on semiclassical arguments. It is well known 
that the number of quantum levels with energy less than E, the counting function N[E), 
is for any classical bounded Hamiltonian Hci{x,p) in one dimension given by (see e.g. 
[32])) 

N{E) = area(E)(l + 0{h)), (29) 
ZTcn 

where 

areaiE) := J dxdp e{E - Hd{x,p)) (30) 

V 

is the phase-space area under the contour Hci{x,p) = E. Obviously, there is a problem if 
this formula is applied to the Hamiltonian ([T]), since the classical motion is not bounded, 
so that a.Tea{E) is infinite. Therefore, Berry and Keating [13] proposed to regularize 
the system by a suitable truncation of phase space in such a way that area(i?) becomes 
finite. 

The regularization proposed by Berry and Keating |13j is to truncate x and p 
by considering the "regularized phase space" V^eg '■= {ix,oo) x (/p, oo) together with 
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the semiclassical condition l^lp = 27rh. This truncation cuts off not only the "small" 
coordinate x < respectively momentum values p < Ip, but it leads for a given "energy" 
E > also to a cut off at the "large" values x = f- respectively P = f- since E = Hc\{x,p) 
holds. Without specifying the behaviour of the classical motion at the end points of the 
trajectories, we follow Berry and Keating and obtain from ( l29l) and f l30l) 

1 



X \ 



'l + 0{h)) 

(31) 



:^ ( In ( A 1 - 1 1 + 1 + . . . . 



N{E) 



2'Kh V \27ch^ 

Setting h = 1 together with a modification of N{E) by adding — | to the right-hand side 
of ( 13T]) which was suggested by Berry and Keating [T3l HI] in order to take into account 
the Maslov index, we arrive at the leading asymptotics of the counting function of the 
nontrivial zeros of the Riemann zeta function (Riemann- von Mangoldt formula) 

Following the argumentation of Berry and Keating [T3|[T^ for the modification of N{E), 
we get for the corresponding Maslov index yU = — |. This seems at first somewhat strange 
since there is no magnetic flux or spinning particle given and, therefore, the Maslov 
index should be an integer number as in the case of "normal" quantum systems like the 
harmonic oscillator. We want to mention that there is actually no rigorous argument for 
the choice of the Maslov index (correction) simply by the fact that so far we have not yet 
imposed any boundary conditions on the operator, and in the corresponding classical 
description there is therefore a lack of jump or scattering condition at the end points of 
the trajectories. The scattering conditions in section [16] (example 116. 2p could provide a 
possible remedy for the above mentioned discrepancy of the Maslov index with respect 
to "normal" systems. Furthermore, there is only one possibility in the classical case for 
the behaviour of the particle at the end point of the trajectory if one wants to preserve 
the constancy of the Hamiltonian for all time: the particle must jump from the point 
{f-,lp) to the point (Z^., ^) in phase space, which corresponds to a kind of ring-system 
(one-dimensional torus with the topology of S^) in the configuration space. 



4. Classical dynamics and quantization of the "squared" Berry-Keating 
operator 

In order to allow some kind of reflection at the end points of the trajectories, we shall 
also consider the classical Hamiltonian 

Hci{x,p) := xV, (33) 
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which is the square of the Berry-Keating Hamiltonian ([T]). Note that fl33p can be derived 
from the Lagrangian 

L{x,x) = \(^iy (34) 
and that Hamilton's equations do not decouple in this case as in ([2j). In fact, one obtains 

x(t) = = 2x^p(t) and p(t) = = -2xp^{t), (35) 

op ox 

and the solutions are 

x{t) = xoe^^o?'"* and p{t) = poe"^"^"^"*. (36) 
If one broadens the phase space to 

'Preg,b ■= {Ix, oo) X ((/p, Oo) U (-/p, -Oo)) (37) 

one now has the possibility to scatter from the end point {f-,lp) of a trajectory of the 
form ( !36|) to the end point (p, —Ip). This corresponds to a reflection on a wall like in 
a one-dimensional billiard system. This is one reason why we rather consider Hd and 
accordingly Hbk as a momentum (operator) and Hd and respectively i^BK energy 
(operator). Further hints to this choice will follow in the sequel. 

Before investigating the "squared" Berry-Keating operator on quantum graphs, we 
would like to consider this operator in the framework of standard quantum mechanics 
restricting ourselves, however, to the positive half-line M> as in the discussion of the 
original Berry-Keating operator in section [2l A formal calculation of H := H^^^ gives 
(setting from now on h = 1): 

HIk := + -) V = - - -■ (38) 

\ \ dx 2/ J dx^ dx 4 ^ ' 

Again as in section [21 i^BK acting on D(]R>) is essentially self-adjoint, and in the 
following we always consider the self-adjoint closure of this operator. It is worthwhile 
to mention that the squared operator (138|) is a special case of the famous Black-Scholes 
operator |33l IM] introduced to determine the pricing of options in financial theory whose 
interesting mathematical properties have been discussed e.g. in [27 1 128 | |29 ] . 

It is easy to see that the functions ipk{x) (A; G M \ {0}) defined in (fTOj) are the only 
eigenf unctions of i^BK ^iii Il^> corresponding to the continuous spectrum A = A;^ > 0. 
Here the eigenvalue A = (respectively k = 0) corresponds to the two generalized 
eigenf unctions 

1-1 1-1 

'^o,i(^) — and ^o,2{x) = x_^_'^ Inx. (39) 

v27r ' v27r 

An eigenvalue A = A;^ > possesses the two linearly independent generalized 

eigenfunctions tpk{x) and tp^ki^)- 

Introducing the (retarded) integral kernel of the time-evolution operator (unitary 

group) 

U{t) := e~"^BK (40) 
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by 

:= (x) =: J ir(x,xo;t)0(xo)dxo, (41) 

R> 

where (j){x) G (M>, dx) is the initial wave function at t = 0, we obtain (cf. [30< p. 27]) 

oo 

K{x,xo;t)= [ Mx)^k{xo)e-"'"'e{t)dk 

I (42) 

1 . (In x — ln xri)^ 

= {AttUxxq) 2e' J^e(t). 

The kernel K satisfies the inhomogeneous time- dependent Schrodinger equation 

- HIk,.^ K{x, xo; t) = i5{x - x^)5{t), (43) 

i.e. it is the retarded Green's function. With (HTl) . the action of U{t) on G (I^>) dx) 
is given by (t > 0) 

oo 

U{t)<p^ (x) = (47rit)"5 / e'i^e-50 {e'^x) dr, (44) 



which expresses the fact that U{t) is a combination of the scaling transformation 
generated by the operator D = x^ (see eq. fl22|) ) and the transformation generated 

generated by T reads 



by the operator T := x^^, since U{t) = e^4e^*^e^'*^. Notice that the transformation 



oo 

e~^4 



(e'*^0) (x) = r / e%+^0 (e-'x) dr, (45) 

(47rit)2 J 

— oo 

and that the operators D and T commute. The resolvent kernel (outgoing Green's 
function) of is given by (see [30, p. 26]) 



G{x,xo;X) :=i J e'^^+''^'K{x,xo;t)dt 



= {AXXO (-A - ie))"5 g-(-A-ie)^|lnx-lnxo|^ 



(46) 



which shows that G has a cut on the positive real axis in the complex A-plane (if ^/z 
is defined with a cut on the negative real axis in the z-plane). With k := \/X > one 
obtains {x,Xq G ]R>) 

G {x,xo;e) = —-^e''\'^^-'^^o\ 
Ik^xxq 

_ iTT J ^/'fc(x)^fc(xo) for x>xq ^^^^ 
k [ ?/'^,(x)V'fc(xo) for X <Xo 
in agreement with the general form of the Green's function of a Sturm-Liouville operator 
(see e.g. [Ml p.ll2]). 
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5. Semiclassical estimate for the eigenvalue counting function of the 
"squared" Berry-Keating operator 

Using again the semiclassical formula (129|) and the truncation of phase space as discussed 
in section [3l we obtain for the counting function in the quadratic case 



where we have included the same Maslov index correction as in section |3l Furthermore, 
we have introduced the "wave number" k, E =: ffk"^, and have used l^lp = 27rh. We 
note that in this case we obtain twice the counting function of the Riemann zeros (for 
which only those with positive imaginary part are counted), since each energy value 
E comes with two values ±k. Notice, that in this case the Riemann zeros are not 
interpreted as "energies" but rather as "momenta" hk respectively "wave numbers" k. 
Formula (HHll agrees with the well-known universal law that N{E) for a bounded system 
in d dimensions grows asymptotically as N{E) = 0(-E'2 ), and thus for a one-dimensional 
system one expects N{E) = 0{y/E) = 0{k), eventually modified by a factor ln(-\/E). 

6. Compact graphs 

We shall present a short overview on compact graphs using the notations of |36] and 



A compact graph F = (V, S, I) is a finite set of vertices V = [vi, . . . , vy) and a finite 
set of edges S = (ei, . . . , e^;). Here we have defined E := \S\ and V := |V| for the total 
number of edges and vertices, respectively. Each vertex f G V is at least connected with 
one element v eV hj some edge e E S, where v = v is allowed. Furthermore, each edge 
e E S connects two vertices v and v in V, again v = v is possible. The topology of the 
graph is given by these relations of the edges and the vertices. Each edge e is assigned 
an interval Jg = [a^, b^] with < < < oo. The set of all intervals is denoted by I. 
We remark that the choice of the starting point and the final point be of the edge e 
is arbitrary and there is no orientation of the graph assumed. We denote two edges as 
adjacent iff they share at least one vertex as endpoint. We need the notion of a path and 
of a periodic orbit of the graph. We slightly differ from the definition in [37j for further 
convenience. A path p{w, z) := {{ei)^^^, w, z) is a set of a finite sequence of edges (ej)"^;^ 
where the points w, z G I denote the starting and final points of the path. Furthermore, 
it is required that 

• the edges Cj and Cj+i are adjacent, 

• the point w must be an element of Ji and z must be an element of /„. 




(48) 
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The case w = z is admissible and corresponds to a closed path. In |37| or [TB] only the 
first item is required for a closed path at which we set Ci = e„+i. We shall call this case 
a closed orbit. Especially, a closed orbit is only characterized by a sequence of edges 
(ci)r=i- ■^o'^ definition of a periodic orbit 7, we shall keep with the usual definition 
as in |37| . then 7 is an equivalence class of closed orbits and can be characterized by 
a representative 7 = (ej)[L]^. The set of all periodic orbits is denoted by ^i. We could 
then equip the graph with a metric structure in an obvious way like in j3^. Especially, 
this would mean that the length of the edge Cj will be U = hi — ai. However, here we 
take another choice for the lengths and the metric structure of the graph. We define the 
length z) of the path p{w, z) := ((cj)^^]^, z) as follows. If n > 2 we denote by yi 
and yn the endpoints of the intervals Ii and /„ corresponding to the shared vertices of 
the edges 61,62 and 6n_i,6„. In particular this means that yi is identical with ai or hi 
and yn is identical with a„ or 6„. Then the length [p(w, z) is defined as 

n-l 



ip{w,z) :-- 
Similarly, if n = 2 {yi = 



lp{w, z) 



In 



In 



yi 



w 



+ In 



hi 



+ 



In 



y) respectively n = 1 we define 



+ 



In 



respectively 



lp{w,z) 



In 



Furthermore, we define in a natural way the length of a periodic orbit 7 

hi^ 



■7 



In 



n 



a,; 



(49) 



(50) 



(51) 



In order to define a metric structure of the graph, we need the notion of connectedness. 
We define w and z as connected iff there exists a path p{w,z) := {{ei)^^-^^,w,z). The 
graph r is connected iff all points of the intervals I are connected. Not necessarily but 
for convenience, we assume in the following that the graph F is connected. The distance 
duj^z of two points w and z on the edges of the graph is defined by 

d^^z '■= niin {^p] P connects w and z} . (52) 

We remark that this choice of the metric for the graph will correspond to a "hyperbolic" 
metric in one dimension. In this case the determinant of the metric tensor g at the 
point X is {detg){x) = The reason for our choice of the metric will be explained in 
sections [7] and [131 

We also need for the interpretation of the trace formula for the Berry-Keating 
operator H^y. in theorem 115.21 the notion of a directed graph in order to interpret the 
right side of (I136p as a sum of periodic orbits. Therefore, we replace the edges by directed 
edges. This doesn't affect the lengths of the edges but has of course an influence on 
topological properties of the graph such as connectedness and on the set ^ of periodic 
orbits. Since in this case there are only such paths p(w, z) := ((6j)[L]^, z) allowed for 
which for all consecutive edges 6j, 6^+1, there exist vertices Vij such that the direction of 
6j is towards Vij and e^+i has the direction away from Vij. Then, the definition for the 
periodic orbits and for connectedness for directed graphs are the same as for undirected 
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graphs. Again, not necessarily but for convenience we always assume that the directed 
graph (in the case for the Berry-Keating operator //bk) is connected. 

7. The Berry-Keating operator on compact quantum graphs 

We define, in accordance with |37] : 

E 

C^i'^)-=®C^KH and 



i=l 
E 



(53) 



n := L'^{T) := ^L'^ {[ai,bi],dx) with < < 6^ < oo. 

i=l 

This means that a function of the Hilbert space "H is represented by an orthogonal sum 
of functions which are defined on the corresponding edges: 

E 

ijjeL\r) iff ^ = 0^i with -^ie L\li,dx). (54) 

i=l 

The first function space in ( l53i) will be a possible operator core for the closed 
Berry-Keating operator with (157|) as domain of definition. Therefore, the self-adjoint 
extensions of Hbk are also with respect to C^(r). The space H is a Hilbert space if we 
equip it with the scalar product 



(^,0) -.= ^2 [ 

•J n ^ 



i)j{xj)(t)j{xj)dxj. (55) 



We then define the Berry-Keating operator on compact graphs (in the following we set 
h=l): 

H^K^P := + . . . , -i(x^ + D^Pe) , (56) 

for ip G C^(r). Since we have a compact graph F, multiplication by x is a bounded 
closable operation. Thus, by perturbation arguments, see e.g. [38| p. 183], we conclude 
that Hbk is closable since the standard momentum operator p = — i^ is closable. 
Furthermore, we note that multiplication by the argument is also a (bounded) bijection 
from 

E 

Dl{r):=^H',[a,,k] (57) 

i=l 

to itself. H^laijbi] is the set of absolutely continuous functions on [aj,6j] which vanish 
at the endpoints of the intervals and with square integrable weak derivatives. Again, 
by perturbation arguments for the momentum operator, we therefore conclude that the 
domain of definition of the closure of Hbk is equal to (157|) . Furthermore, by similar 
arguments the adjoint operator of (i^BK, -Do(F)) is given by {Hbk, H^(r)) in which 

E 

H\r) :=^H'[a„k] (58) 

i=l 
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is the set of absolutely continuous functions on the intervals of the graph F, cf. |38[ 
p. 100] possessing square integrable weak derivatives. 

We mention at this point that the projections of the spaces DKT) and H^iT) on 
the intervals of the graph F coincide with the corresponding Sobolev spaces, see e.g. [39j. 
The operator (i^BK, -Do(F)) is symmetric and it is possible to show that the deficiency 
indices are {E,E), compare e.g. p. 142]. By a proper Sobolev embedding theorem 
and the compactness of the graph F it follows that the differential operator on -Do(F) 
possesses a compact resolvent, see also [40j. Thus, by the compact resolvent theorem, 
cf. 1^ p. 245], and the relatively compact perturbation theorem (cf. [1T| p. 113]) the 
operator (l56l) possesses a purely discrete spectrum. 



8. Classification of the self-adjoint extensions of the Berry-Keating operator 

In order to characterize the self-adjoint extensions, we follow the ideas of [24, P- 138] 
and [36j . see also |42) for a comprehensive discussion. Therefore, we define the complex 
symplectic form on H\T) x H\T) (cf. [Ml P- 138]): 

[0,^]i:=(0,i/+KV^>^,(^)-(ff+K0>V^>i2(r) for <f),^ e H\r). (59) 
We call a subspace X [-, -j-^-symmetric, iff [0, '?/'];^ = for all G X. Due to the 
von Neumann extension theory (see e.g. [21]) the self-adjoint extensions are exactly the 
maximal [■, -jj^-symmetric subspaces of if^(F). We follow the approach by Kostrykin and 
Schrader [36] to classify these extensions. By a proper Sobolev embedding theorem we 
can define the boundary value bv (an element of C^'^) of & H^{r): 

^i,.:=iM(^i),...,MO'E),Mbi),...,MbE)f for e H\T). (60) 
For convenience, we also define: 



tExE 
~1e>cE 



'± - I n T ) ' ^(-f) •= I n I. I ' (61) 




with (no summation over i) 

O'ij '■= Sijai and bij := 6ijbi for < i, j < E (62) 

and 

J:=l ° ^^^^ 1, U:=^\ 1. (63) 




-Iexe 

By a simple calculation we obtain the identity: 

f/(i/±)f/+ = J. (64) 
Thus, we obtain for f l59|) by an integration by parts using the unitarity of U : 

= ($bv,i/±^(a6)^bv)c2i5 

:UDl,^'^>^.,JUDl,^^^^ for all V^,0Gif^(F). 
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To define the square root of D^ab) we liave used the usual definition of a positive operator 
( |39[ p. 196]), which in this case simply means to take the square root of the (diagonal) 
entries in D(^ab)- Note that 

defines a nondegenerate complex symplectic form on x C^^. We call a subspace C 
of C^"^ a Lagrangian subspace iff 

• a,b E C then u{a,b) = 0. 

• Whenever for a subspace C D C the first property holds, (67) 
it follows C = C 

For the Lagrangian subspaces of C^^ we apply the result of [3^. A subspace C is 
Lagrangian iff there exist two matrices A,B E Mat(i? x E, C) with: 

AB+ = BA+ and 

(68) 

rank(A, B) = E. 

We then have, 

G C^^; := ^ j and A0i + S02 = o| . (69) 

In ( l68l) the matrix (A, B) is formed of the columns of A and 5, and we have introduced 
two maps 

(.)^:C2^^C^ for l<i<2 (70) 



(1,0)0 ifz = l, 
(0,1)0 if 2 = 2. 



by 



Furthermore, as mentioned in [l3], these matrices are not uniquely defined. Two sets of 
matrices A, B and A, B define the same Lagrangian subspace iff there exists an invertible 
matrix C with 

A = CA and B = CB. (72) 

1 

Since UD'^^^-^ is a bijection from C^^ onto itself and by (159|) and ( 165|) we infer, with 
the same arguments as in [36], that there is a one-to-one correspondence between the 
self-adjoint extensions of {H-q]^, Dq{T)) and the Lagrangian subspaces of C^^. We thus 
have proved the following proposition. 

Proposition 8.1 Each domain of definition of such a self-adjoint extension is exactly 
the preimage with respect to / fgOj) of a subspace 

^ = {^bv G C^^; A [udI,^^>^:)^ + B (t/I^(L)^bv)^ = O} , (73) 
where A and B fulfil [6B\] . The converse is also true. 



Because of proposition 18. ![ we denote in the following the self-adjoint extensions of 
{H^^,Dl{V))hj {H^^-A,B). 
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9. Determination of the eigenvalues of H^k 

All possible eigenf unctions tpk to an eigenvalue k of Hbk are of the form 

M^) = f al^e''='°^ . . . , aE^e'^'^A . (74) 



We denote the column vector fl60l) corresponding to ipk by '^bv,k- In order to apply (173|) 
for determining the eigenvalues k and the corresponding eigenvectors ipk, we calculate 
UDl^^^bv,k using dMD. 

1 / gifclna \ 



e 



1 / jgifcliK^ gifclnb 
1 / jgi'^lfi'^ _|_ gifclnb 



a (75) 



a. 



ik In am 

(76) 



For convenience, we have used the notations: 

- 5mne'^'"'™ for 1 < m, n < E 
and 

CK := (ai, . . . , a^, ai, . . . , a^;)"^. (77) 
Therefore, we get: 

) \ , (78) 

[UDfab)'^iro,k)^ = - (e^'^"" + ie^'^"') "2, 

where we have used (•)j defined in (1711) . 

Taking into account that tXi = a2 ='■ ck, we obtain for the expression in (173|) 

= (^A(ie''=^''" + e''^'°'') -5(e^^''^'^ + ie^^'''^))Q: (79) 

= (i (A + 15) e^'^^"" + (^ - iB) e'^'"^'') a. 

Kostrykin and Schrader have shown, see [36j, that A ± iB are invertible under the 
assumption ( 168|) . With the notation 

C(A;) :=e'*^^""a (80) 

and, because {A — 15) and (A + iB)^^ commute (see [44j), we get: 



^(l-i^e--i^)C(A;)=0, 
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with e''^''^a similarly defined as in (1761) . Due to the similarity of (IHTl) with the secular 
equation for the common Laplace operator on compact graphs (see |l3]), we denote: 

5(A,5):=i4— ^ and T (a, 6; fc) := e''^^" ^ (82) 

It follows with exactly the same arguments as in [36j that S{A, B) is unitary, and we 
shall call it also the 5'-matrix of the quantum graph. The unitarity of T (a, b; k) iff 
/c G M is obvious. Since, for all A; G C, C(/i;) = iff ck = 0, we have proved the following 
proposition establishing the secular equation J-'{k) = 0. 

Proposition 9.1 A;„ G M an eigenvalue of {H-q]^\ B) iff 

7{kn) := dei{lExE - S{A, B)T (a, b; k^)) = 0. (83) 

Furthermore, the multiplicity of the eigenvalue one of S{A, B)T {a, b; kn) coincides with 
the multiplicity of the eigenvalue kn of {H^i^; A, B) . 

We remark that the geometric multiplicities and the algebraic multiplicities of 
S{A, B)T (a, 6; k) coincide since this matrix is diagonalizable. 

In contrast to the S-matrix of the generic negative Laplacian —A on graphs, the 
S-matrix S{A, B) is always independent of k (the S'-matrix of —A is independent of 
the wave number k iff = S [see |13]]). But we remark that the independence of 
the S'-matrix on the eigenvalue will also occur when we replace Hbk by the standard 
momentum operator (with x G M) 

P ■■= -i^. (84) 
ax 

The calculations are quite analogous (see e.g. [45] for a detailed discussion of this). In 
fact every self-adjoint extension of p can be characterized by the same matrices A and 
B as in (168|) and we would get the same S'-matrix S{A, B). The only difference in the 
secular equation between the operators i^sK and p then is the form of the second matrix 
in (l82l) which in the case of p is given by 

T(a,6;fc) := e^'=(''-"). (85) 

This is one reason why we rather relate i^BK with a momentum operator than an energy 
operator as indicated in section [2J Because of the occurrence of the logarithm in T we 
endow the quantum graph with a metric structure as proposed in section [6l We also 
mention that proposition 19.11 is valid for all G C, in particular for = in contrast to 
the corresponding proposition 112.11 for H^j^. However, the analogy between p^ = —A 
and -f^BK ^'^ obvious. 



10. The "squared" Berry-Keating operator 

Our Hilbert space will be "H = L'^{T), see (15^ and ( 155|) . Again, we seek self-adjoint 
extensions of fl38l) with respect to C^(r). In order to obtain a self-adjoint operator, the 
task is to specify an appropriate domain D[H^^) for this operator with 

C^iT) C D{H'^^). (86) 
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One simple possibility is to define H"^^ as the "squared" Berry-Keating operator, which 
means: 

Hl^^l^ := HnK (Hsk^P) , 

^ e D{Hl^) := {0 G DiH^K); HBK<f> G D(i/BK)}. 

It follows immediately that -^bk self-adjoint if i^BK is self-adjoint using Friedrichs' 
extension theorem [Ml p. 180]. But in fact there are many possible self-adjoint extensions 
which cannot be realized in such a way. We will give simple examples in section [161 We 
can generalize these constructions to consider non-self-adjoint but closed realizations of 
Hbk and then form 

i/+K^BK or i^BK^^^K- (88) 

This is an idea quite analogous to the concept of supersymmetry, see |46] and |47] (the 
technique of factorization was already introduced by Schrodinger [l8] and reviewed in 
[49] ) . In |17] this technique has been used but isn't explicitly mentioned. However, 
only a certain kind of self-adjoint extension can be attained in such a way. In [17] these 
are exactly the self-adjoint extensions which correspond to fc-independent S-matrices 
corresponding to these extensions. This relation between the S-matrices and the self- 
adjoint extensions of the negative Laplace operator —A on metric graphs is explained 
in [43]. 

We would like to give an overview of the starting point of our considerations from 
a mathematical point of view. The proofs of these statements are similar as in section 
[7] using the same references as there. Therefore, we only summarize the results. First, 
the operator -f^BK acting on C^{r) or 

E 

^o(r) := H^K bi] with < a, < 6, < oo (89) 
1=1 

is symmetric. H^[ai,bi] is the set of absolutely continuous functions which possess 
absolutely continuous derivatives on [aj,6j], square integrable weak second derivatives 
and which together with their first derivatives vanish at the endpoints of the intervals. 
Furthermore, i^^K acting on Dq^T) is closed and the adjoint operator of {H^^jDqIT)) 
is {H^^,H\r)). Here 

E 

H^{r) ■= H^[a„ hi] with < < 6^ < oo (90) 

1=1 

is the space of functions being absolutely continuous on the intervals of the graph 
r possessing absolutely continuous derivatives and weak square integrable second 
derivatives. The deficiency indices are {2E,2E), thus (H^^jD^lT)) possesses infinitely 
many self-adjoint extensions. Again, as for Hbk the spectrum of every self-adjoint 
extension is purely discrete and as in section [7] the projections of the spaces Dq^T) and 
i7^(r) on the intervals of the graph F coincide with the corresponding Sobolev spaces 
(see again e.g. [39j). We shall follow a general approach to find all these self-adjoint 
extensions, quite analogous as in section [8] and based on [36j. 
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11. Classification of the self-adjoint extensions of the "squared" 
Berry-Keating operator 

First, we define \[^bv as in f l60|) for G -f^^(r) and additionally 

^bv := i^'Mi), • • • , ^'eM, -tp'iibi),. . . , -tp'E{bE)f 
for ^ e H^iT), 

in which ■^|)^ is the derivative of ijji on the interval Jj. Similarly as in (l59l) . we define a 
symplectic form on H'^{T) x i^^(r) 

2 +„/,\ / Tj2 +, „/,\ r„„ , „;. ^ u2 



(91) 



[0, V'ls := (0, ^BK>;^,^^^ - (^BK>, for 0, ^ G H\T). (92) 

With the same arguments as for Hbk in section [7] the task is to find all maximal [■, ■]2- 
symmetric subspaces of H'^{T) in order to find all self-adjoint extensions to {H^-^, Dq^T)). 
We shall adapt the definition of J in (163|) by 

J := I ° ^'l^'"^ I (93) 



and (see (I6T 



We obtain for 4>,ip E H^{T) using partial integration and the fact that J and -D(ab) 
commute: 

E 

1=1 

Mbv,j5?„,)[0]bv)^^^ ^^^^ 



(96) 



^(a6)[^]bv, ^^(a6)[0]bv, 

Taking the scalar product in the definition of a;(-, ■) in (1^^ with respect to , we infer 
as in section [8] the following proposition. 

Proposition 11.1 The self-adjoint extensions of (iJgj^, Dg(r)) are exactly the 
preimages of 

C = |[0]bv G C^^; A (5(„5)[0]bv)^ + 5 (5(„5)[0]bv)^ = O} 
= {[0]bv e C^^; AD(„5)$bv + 5D(a5)$bv = 0} 

with respect to [-Jbv [94\ )- 

In fl96|) we have used (•)j defined in fITT]) . The matrices A and i? are now elements of 
Mat(2£' X 2E,C) with the adopted conditions 

AB^ = BA+ and 

(97) 

rank(A, B) = 2E. 

Again, as in section |8] and because of proposition 111.11 we denote the self-adjoint 
extensions of {H^j^, D^{T)) by {H^^;A,B). 
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We want to solve the eigenvalue problem 

H^K^ = X^. (98) 

To tackle this problem, it will be convenient to consider the wave number k defined by 
A2 := k. It is a trivial observation that ±k correspond to the same eigenvalue A. This 
fact will be revealed in the symmetry of the secular equation for the wave number. 

In addition to 098 p the eigenvector ip must be in the domain of definition of the 
operator. However, the general form of the eigenvector to an eigenvalue A = fc^ 7^ is 



(„ „iA;lna: i n „— ifclna;\ /„ „ifclna; i n ^—ik\nx\\ /'nn^ 

(aie + Pie j , . . . , (a^e + p^e -(99) 



We can proceed as in j36]. Therefore, we compute '^bv,k and ^E'^^ using the definitions 
in (ETl) and (176]). 



gifclna g— ifc In a 



^bv,k — -0(a6) \ gifelnb g-ifclnb 




1 3 / pifelna — ifelna 

= I -2^(al) _gifcln. _g_ifcln. 1 (10°) 




gifc In a 
+i^-^(al) [ _gifcln6 

In order to be in the domain of definition of a self-adjoint realization, '^bv,k and ^^'^^ must 
be in some C of ( p6|) defined by the two matrices A and B. We make the identification 
(cf. 03]), 

pifclna p— ifc In a 

Xik;a,b):= \ |, (101) 

pifclna p— ifc In a 

Yik;a,b):= | ) and (102) 

pifclna p— ifc In a \ 

Y'ik;a,b):= \ \. (103) 

Thus, we conclude, with the definition for the bold symbols in accordance with f l62|) : 

= AD(,5)$bv + 5^(a5)$'bv (104) 




((^^L) - \BD;i)I^)x{k, a, b) + \kBD-i^Y'{k- a, 6)) ^ 
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At this point we make two observations: Since is self-adjoint, we conclude 

ADl,^{BD^l)y = AB^. (106) 

Since -D^^^^ and D^^^^^ are invertible and diagonal, it is easy to show that 

rank(ADf„,), = rank(A, B) = 2E. (107) 

Therefore, we define 

AdI,)=--^ and BD;i^=:B' (108) 

and observe that A' and B' also fulfil the conditions (p7|) . Therefore, we can apply a 
theorem of Kuchment ^40j. It states that two matrices A' and B' fulfil (l97j) iff there 
exists an invertible matrix C with: 

A! = CP^,,B' + CPi,j,,L'Pi^j,, and B' = CPi,j,,. (109) 

In ( llOOp we have defined PkerB' as the projector onto the kernel of B' and -Pj^i-b' 
corresponding orthogonal projector. The matrix V is self-adjoint and can be defined by 
(see [40J), 

^':=(5lanB'+)"'^'nirB'- (HO) 

Hence, we can proceed in the calculation (I104p by multiplying f llOSp from the left-hand 
side by C^^ 

+ ifcP45,r(A:;a,b)j f ^ 

Since the projectors PkcrB' and Pj^j-b' mutually orthogonal, we infer from fll09p . the 
definition of C in (p6|) and with the first line of (llOSp . that 

PkerB'X(fc;a,6) 1^ ^ j =0. (112) 

In ( imp we insert between the matrices I± and X{k; a, h) the unit matrix 1 = 
PkerB' + Pkiij' and apply ([m]) 

= ((PkcriJ' + Pkiri.' - \l±)PLB)x{k\ ft) 

(113) 

+ \kPt^,Y'{k-aM \ \ ^ 



We realize that L' — \l± is also self-adjoint. Thus, we define 

1 
2 



L" := P^^, ( L' - ] P^t^, (114) 
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and make a re-definition: 

A" ■.= P^,,B' + L" and 5" := P^irS'- (115) 

It is obvious that the matrices A" and B" fulfil the conditions ( 195|) . Hence as in section 
|9] respectively [36], we infer that A" ± ikB" is invertible and conclude with quite the 
same calculation as in |13] 

= {A" + ikB") [1 - S"{A, B; k)T{a, b, k)] I „ Jin5 ]i Z \ - (US) 



e-^'^^'"" / \ f3 



Here we have used the definitions 



A" — \hR" 

S"{A^ B- k) := S{A\ B"- k) := and 



n e 
na,h-k) := I ^^^^^ 



ik In 



:ii7) 



e"'^"^ 

The first and the third matrix in the product of (I116p are invertible for all 
A; G C \ (±icr(L")) in which cr(L") denotes the spectrum of L" . For a detailed discussion 
of this, see |37| and [44j. Thus, we have proved the following proposition. 

Proposition 12.1 k'^ with k e C \ (±icr(L") U {0}) is an eigenvalue of {H'^x'^ A, B) iff 

F{k) := det il2Ex2E - S"iA, B; k)T{a, b, k)) = 0. (118) 

Furthermore, as in section [3 the multiplicity of the eigenvalue X = k"^ coincides 
with the multiplicity of the eigenvalue one of S"{A, B; k)T{a,b; k) for every k G 
C\ (±ia(L") U {0}). 

We remark that the restriction on k concerns only non-positive eigenvalues A = fc^ of 
13. The eigenvalue zero 

For the eigenvalue A = 0, which is equivalent to the case k = 0, the eigenfunctions are 
of the form 

V'o(a^) = ( + A^lnx, . . .,aE-]= + /^^^Inx] . (119) 



With a similar calculation as for the case k ^ Q one obtains the equation 

/ 1^.^ Ina V f ^E.E \\(A=Q (120) 

which is necessary and sufficient for A = to be an eigenvalue of {H^y^;A,B). The 
matrices A" and B" are the same as in f lllSp . Then we can proceed as in [37] and get 
the following proposition. 
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Proposition 13.1 X = = is an eigenvalue of (-f^BK' ^) '''ff f'^'^ ^'^^ value k' ^ 
and then for every k' 

Fo{k') := det(l - S"{A, B; k')C{a, b, k')) = (121) 

is fulfilled with S"{A, B; k') as in jjllT]) and 

C{a, b; k') : = 







2f 



2 -V +ln 



k' ' \aE) 







2f 



2f 



-In 



v 







2f 







In 



(122) 



Furthermore, the multiplicity of the eigenvalue A = coincides with the multiplicity of 
the eigenvalue one of S"{A, B; k')C{a, b; k') for every real k' ^ 0. 

Thus, in general there is a difference between the spectral multiplicity of the eigenvalue 
one of S"{A, B;0)T{a,b,0), which we denote by N, and the eigenvalue one of 
S"{A, B- k')C{a, b, k') with k' 7^ 0, see |37], HZ] and p]. 

In order to relate the self-adjoint extensions of {H'^y^, Dq{T)) with the self-adjoint 
extensions of the Laplacian —A one has to adjust the lengths as before. However, in 
oder to attain the same spectrum, except for the case k ^ C \ {±ia{L") \ {0}), one has 
to transform the matrices A and B into A" and B" as in ( llOSp and (IllSp . Then the 
spectrum of the negative Laplacian characterized by A" and B" with the previous choice 
of the lengths will coincide with the spectrum of characterized by A and B in (p6|) . 
Especially, the functions F{k) and Fo{k) in flllSp and (I12ip . respectively, will coincide 
with the corresponding functions for —A, see e.g. |37j and |43j. 

We remark that the transformation of the matrices A — )• A" and B — )• B" and 
vice versa (actually we consider below the converse direction) cannot be achieved by 
perturbing the negative Laplacian by a magnetic flux which corresponds to an operator 
acting on the edges as 

2 



dXj 



iAj{xj^ 



for l<i<E. 



(123) 



Kostrykin and Schrader have shown in [51] that this operator is related to the negative 
Laplacian —A by a unitary transformation of the corresponding S-matrices. This means 
that the Laplacian perturbed by a magnetic flux can also be characterized by two 
matrices A and B obeying f l971) . But with a local gauge transformation this system 
can be transformed to a quantum graph system with the pure Laplacian which is now 
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characterized by two new matrices A and B. These new matrices are obtained by the 
old ones by 



where [/ is a diagonal unitary matrix. If we calculate the S-matrices for these systems 
we obtain 



In particular this means that the S-matrix is /c-independent iff the original S-matrix 
is fc-independent. By a result of [52] we conclude that in the sense of (I109P (see |40| ) 

the corresponding matrix Pi^j-B^PkerB ^kcrB'^'^kerB ^^^o- This feature is 

obviously not given by the transformation A, B to A", B" especially in (I114p taking 
into account that the transformation A, B to A', B' in fllOSp and the corresponding 
transformation L to L' possess this feature. 

Kostrykin and Schrader have shown in |44) that the negative Laplacian —A 
possesses time-reversal symmetry iff S'^ = S. Obviously, the transformation (I125P 
doesn't maintain this symmetry in general. 

In both cases {Hbk and H'^j^) we get the same length U for the edge Ci of the 
quantum graph for the corresponding momentum operator or kinetic energy operator. 
Thus, we choose U for the lengths of the graph and endow it with a metric structure as 
in section [61 

14. Connection between Hbk and H'^j^ 

We also want to reveal the link between Hbk and H]^j^, the last one considered as 
the "squared" Berry-Keating operator as in ( IHTl) . Therefore, we are starting from 
[Hbk] a, B) with corresponding S-matrix S{A, B) and then calculate [H'^^; A, B) with 
corresponding S-matrix S{A, B). Using the definitions (|73l) we obtain from ( !87l) the two 
(necessary and sufficient) equations, setting (p := Hbk^^, 



A = AU and B = BU 



(124) 



S{A,B;k) = US{A,B;k)U+. 



(125) 




(126) 



Making the definitions 



7/^(0) = (V^i(ai), ... ,V'i7;(a£;))^ and ipib) = {ipi{bi), ... ,'ipEibE)f (127) 



the equations ( I126p can be transformed into (see definition (162|) ) 




(128) 




3 



The Berry-Keating operator on (M>,da;) and on compact graphs 25 
This is equivalent to the equation 



1 S(A.B) IW.-O. 



(129) 



Comparing f ll29p with (p6|) . we infer for the matrices (one possible choice) A and B 

~ fo \ 

. > ^ = /^x^ • (130 

It is a simple calculation that indeed the matrices A and B fulfil (p7|) . In order to 
calculate S{A, B) we make two observations. First, we infer from fll06p and (11301) that 
L' = in the decomposition (I109p . Furthermore, we notice that 

~,± { ( S{A,B)e, \ ■ , J, 
ker B = span < | | ; 2 = 1 E 



:i3ii 



ker5' = span I '^^^jf j ; ^ = 1,...,E 

where span denotes the linear span of the corresponding vectors, and the vectors Cj span 
the space C^. Therefore, by (I114p we conclude 

L" = L' = 0. (132) 

Thus, by (I117p and after an easy (but lengthy) calculation we obtain for the S-matrix 
of {Hl^;A,B) 

Therefore, we have proved the following proposition. 

Proposition 14.1 {H'^j^; A, B) is the "squared" operator of some {Hbk', A, B) in the 
sense of l[8l\) iff the corresponding S-matrices fulfil U33\) . 

We remark that a similar relation holds for the usual Laplace operator —A and the 
momentum operator p on compact graphs. 



15. Trace formulae and Weyl's law 

We are now in the position to give explicit expressions for the behaviour of the eigenvalue 
counting functions for large eigenvalues and give trace formulae for the Berry-Keating 
operator and the "squared" Berry-Keating operator on compact graphs. These results 
are immediate consequences of sections [TJ [10] and the results in [37]. The proofs of the 
claims for the Berry-Keating operator are quite analogous to [37] and, therefore, we only 
give a short outline of some steps of the proof. Since the trace formulae differ in some 
details we formulate these formulae in one theorem and one corollary. First of all we 
introduce an appropriate space of test functions as in |37] . 
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Definition 15.1 For each r > the space consists of all functions /i : C — )■ C 
satisfying the following conditions: 

• h is even, i.e., h{k) = h{—k). 

• For each h & Hr there exists 6 > such that h is analytic in the strip Mr+s := {A; G 
C; I ImA;| < r + 6}. 

• For each h G Hj. there exists r] > such that h{k) = O ^(^1^1^1)1+?^ on Mr-^.s, 
/c — 7- oo. 

We denote by kn the "energies" respectively "wave numbers" of Hbk respectively -f^BK 
and by the corresponding multiplicities which are identical with the order of the 
corresponding zeros fc„ of J-" in (153]) for n G Nq respectively zeros fc„ 7^ of F in flllSp 
for n G N. n = corresponds to the "energy" respectively "wave number" zero, and the 
energies respectively the nonnegative wave numbers are ordered with respect to their 
absolute value \kn\ in increasing order. However, the (finitely many) imaginary wave 
numbers are omitted. Furthermore, we denote the self-adjoint realizations characterized 
by ( 1731) respectively ( 196|) by {Hbk',A,B) respectively {H'^j^; A, B). Notice that in the 
first case A,B e Ma.t{E x E, C) whereas in the second case A,B e Mat(2£' x 2E, C). In 
addition we denote by Imin the minimal length of the graph with respect to the definition 
of section [6l The minimal positive eigenvalue of L" in fllMp is denoted by X"^i^ and the 
unique minimum of the function 

1 2 / K \ 

i{K) := - \n{2E) + - artanh ( -^jj^ 1 (134) 

mill / 



by a. For convenience, we denote the total length of the graph by 

E 

£:=^[i. (135) 

i=l 

Furthermore, by a hat ■ we denote the Fourier transform (see (ITTl) ) and ■ * ■ denotes the 
convolution of two functions in the distributional sense, see e.g. [24j. For convenience, 
we assume that the graph F is local with respect to the S-matrix S"{A, B; k) which 
means that the scattering between two endpoints is only allowed for adjacent edge ends, 
see [ 44| for a precise definition. This has the effect that in the trace formula the periodic 
orbits are with respect to the classical topology as explained in section [6l Otherwise 
we must interpret the periodic orbits with respect to the topology induced by the S- 
matrices which will differ from the one in section [6] and must then be interpreted as a 
quantum mechanical topology. However, in |31] it was shown that there exists always 
at least one graph, for which the S-matrix is local. In order to interpret the right side 
of fll36p for Hbk as a sum of periodic orbits, we replace the edges by directed edges 
as mentioned in section [6l Furthermore, we assume that the S-matrix S{A,B) is local 
with respect to the directed edges. This means that S{A, B)ij = if Cj and cj share no 
vertex Vij for which cj has the direction towards Vij and Cj has the direction away from 
Vij. Again, it is always possible to find such a graph. We get the following theorem for 
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Theorem 15.2 (Trace formula for //bk) Let T be a compact metric graph and 
(i^BK; A, B) with the above assumptions be given. Let h G with any r > 0. Then the 
following trace formula holds (where h denotes the Fourier transform of h defined as in 
eq. ni\) and Qn the multiplicity of the eigenvalue kn) 

oo 

^gnh{kn)= £/i(0) + 2^Re(^^)/i(^). (136) 

n=0 7e'P 

The amphtude functions are constructed from the S-matrix elements with respect to 
the periodic orbits 7, see |37t [T6] for a precise definition of this construction. The proof 
of this theorem is quite analogous to |37]. Since we have no fc-dependence of 5(^4, 5) 
in (182|) , we can omit the requirement of the minimal length in contrast to the following 
corollary 115.31 This also leads to the simple product of the real part of the amplitude 
functions and the Fourier transform of h in the identity (11360 . Furthermore, since 
the secular equation (18T|) respectively (183|) holds also for the eigenvalue zero of i^BK? the 
term — does not appear in f ll36p in contrast to f ll37p for H^^, where denotes 
the multiplicity of the (possible) zero fco = of F{k). For -f^BK S^^ following 
trace formula. 

Theorem 15.3 (Trace formula for H^^) Let T be a compact metric graph and 
{H'^j^]A,B) with the above assumptions be given. Let the condition [mm > K^) 
fulfilled and let h G Hr with r > a. Then the following trace formula holds 

f2 gn Kk^) = & MO) + (go - \N)hiO) - i- hik) ^^^'^S"{A,B;k) 

+ J2[Ch*A,)H,) + Ch*A,)H,) . 

Again, the amplitude functions A^ are constructed from the S-matrix elements with 
respect to the periodic orbits 7 and qq denotes the multiplicity of the eigenvalue one of 
S"{A, B- k')C{a, b, k') for any A;' G M \ {0} (see section [13]). 

Since we have previously seen that the spectrum of {H'^j^]A,B) coincides with 
some self-adjoint realization of —A on the graph by adapting the lengths and with the 
results in ^7\, we get Weyl's law: 

Theorem 15.4 (Weyl's law for H^^) Given the eigenvalues of some {H'^j^]A,B) 
in increasing order denoted by = k^. Then for the counting function N{X) : = 
7^ {n; k"^ < A} the following asymptotic law holds 

N{\) ~ -V\ for A ^ 00. (138) 

The same asymptotic law holds for {Hbk', A, B) replacing A by on the left-hand side 
and replacing \/~\ by k at the right-hand side on the equation, i.e. we have 

Theorem 15.5 (Weyl's law for i^Bx) Given the positive eigenvalues of some 
{H-BK',A,B) in increasing order denoted by A„ = A;„. Then for the counting function 
N{k) := # {n; kn < k} the following asymptotic law holds 

N{k) ~ -k for k^oo. (139) 
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The theorem 115.51 can be proved by a suitable Karamata-Tauberian theorem (see e.g. 
|53] and [3^) or alternatively by applying theorem 115.41 to the "squared" operator H]^j^ 
of Hbk in the sense of ( IHTl) (the spectrum of the eigenvalues of Hbk and the wave 
numbers of i?^^ coincides then and therefore the corresponding eigenvalue respectively 
wave number counting functions are the same). Comparing the theorems 115.41 and 115.51 
with the asymptotics of the counting function for the nontrivial Riemann zeros (132|) . we 
therefore can conclude: 

Theorem 15.6 (No-go theorem) Neither Hbk nor -f^BK yi'dds as eigenvalues the 
nontrivial Riemann zeros if these are self-adjoint realizations on any compact graph. 



16. Simple examples 

We shall give a simple example for a wave packet and its time-evolution with respect 
to the Berry-Keating operator in = L^(R>,da;) discussed in section O Furthermore, 
we give an example for a realization of i^BK and -^bk the simplest construction of a 
graph which consists of a single edge. Finally, we present some trace formulae for the 
presented examples. 

Example 16.1 For ?/'(x,0) = (j){x) in [2B\) we define (x E M>j 

a 1 
(j)(x) := with a = . (140) 

(With this choice for a it holds = 1.) From ^2^] we obtain 

_ t 
cte 2 

i,{x,t) = {U{t)ct)){x)= with tGM. (141) 

Thus, we get the large-t asymptotics 

01 _t 

■?/' ~ — e 2 for t — )■ oo. (142) 
On the other hand, with 

oo 

KBK{x,Xo;t) = j iJk{x)i^kixo)e-''''dk, (143) 

— oo 

Ijm^ and we get 

oo 

ipix,t)= I A{k)i)k{x)e-'^'Ak. (144) 



A direct calculation using and the integral representation of C,[s) as a Mellin 

transform (see \5^p. 20]) yields 



A{k) = ^ (l - V22^') r Q - ik^ C - . (145) 
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With (see \5^p. 13]) 

r 



27ie for k e 



\k\ — )■ oo, 



(146) 



we get for the large-k asymptotics of \A{k)\'^ 

\A{k)\'^ ~ fs -2V2cos(A;ln2)) e""!^'! 



(147) 



for k G 



\k\ — )■ oo, 



which gives a sufficient condition for A G L^(]R, d/c). // we consider the continuous 
representation (TJ^ of ilj{x,t), we see that tp{x,t) gets no contribution from the wave 
packet A{k) exactly at the wave numbers k corresponding to the conjectured nontrivial 
Riemann zeros. This is reminiscent to the absorption spectrum interpretation of the 
nontrivial Riemann zeros by Connes [771 [7^ . but of course reveals no insight to the 
position of the nontrivial Riemann zeros. 

Example 16.2 For a single edge I = [a,b] (one-dimensional quantum billiard) the 
matrices A and B are arbitrary numbers fulfilling (E^j- The equations ^73) and 
lead then with 



S{A,B) =: e- 



-2nic 



to 



S{A,B 



(148) 



(149) 



27ric 



ip{b) with CG [0,1). 



The eigenvalue spectrum is given by 

271 



kn 



n + c) with c G [0, 1) and n G Z. 



(150) 



In^ 

a 

We now want to calculate as defined in l[8l\) with (TJ^, in particular the S- 

matrix and then compare it with the results in section [7^} In order to distinguish 
the characterizing matrices, we denote these with the subscript -Hbk ^'^^ '^^Ik' ^i'^^'^! 
we derive the transformation of Ah-^^^Bh^^ into Af^2^,Bjj2^ for the corresponding 
operators related by l[8l\) . We get the additional condition 



e 



-27ric„;./ 



ij'{b) with c G [0, 1). 



(151) 
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The two conditions are equivalent to 



-1 i^r 



2g-27ric \ / 



(152) 



Therefore, we define 

f _1 (a\h p-2nic \ / . 

:= ^ « an(^ 5^2 := i ^. (153) 



I ^BK • I ^ (^^y. 



e 



and recognize that indeed 

zs fulfilled. By a comparison of U5S\) with (E^, we infer that ^h^^ '^'^d Bjj2^ are two 
possible matrices to characterize in the sense of (E^. For S"{Ajj2^, Bjj2^] k), we 
get 



V H^^, H^^, J 1 ^2n,c \ S{A,B)+ I ^ ' 



in complete agreement with U33\) and for the secular equation U18\) 

= (^e<'=^^(^)+H _ l) (e<'=^<^)-H - l) . (156) 



J! 



This leads to the "wave numbers 

kn = — ^(n±c), n G Z (157) 

a 

with c as in U50\) . Obviously, with U57\) and U50\) Weyl's law is fulfilled even for 
small n. Alternatively, since these are from the classical point of view integrable systems 
we can perform an EBK- quantization for i^BK o-nd [331 In this semiclassical 
quantization rule the spectrum consists of energies En (for convenience we use now the 
same letter En for kn respectively A„, as in the sectionslE respectively^^ for which (h = 1) 

In{En) = (n + ^) with n>0 (158) 



4 

is fulfilled. Therein fin denotes the so-called Maslov index and 

In{En) = ^ [pdx (159) 



2vr „ 

7n 
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is the classical action of a periodic orbit 'jn which is a subset of the hypersurface H^i = 
respectively Hd = E^. For H(.\ in with the ring system structure mentioned in section 
0, we get from ^EE) 

En = ^(n+^). (160) 



In (^) V 4 

and for Hd in (also with the ring system structure) 

y/E:=K = r^(n + ^]. (161) 



In (!) ^ 4 

A comparison of fll60p with f llSOp yields for the Maslov indices fin = 4c for Hd- For Hbk 
we get two Maslov indices, fin = 4:C for n = 0, 2, 4, . . . and /i„ = —4c for n = 1,3,5,.... 
Since a Maslov index is at most defined modulo 4 and because of c G [0, 1), the above 
second Maslov indices fin = —4:0 correspond to the Maslov indices /i^ = 4(1 — c) for 
n = 1, 3, . . .. We stress that the EBK-quantization for Hd with a classical "hard wall" 
boundary condition yields 

E~n = kn = -^(n + ^), (162) 



In (!) ^ 4 

which differs from (116ip by a factor 2. We mention that the S-matrix elements of 
S"{Ajj2 ,Bjj2 ; fc) for Dirichlet (D), Neumann (N) or Robin (R) boundary conditions 



are given by 



S" { Arr2 , Btt2 I k 

\ ^bk' ^bk' 



—6ij for D, 

6ij for A^, 



Pj - ik (163) 
-Oij \ ., for R, 
Pj + ik 

with Pj e M and i,j G {1, 2} , 

which obviously differ from fllSSp and, therefore, cannot originate from a "squared" 
Berry-Keating operator. If we impose Dirichlet boundary conditions at both interval 
ends, we get for F{k) in f lllSp 

F{k) = l-e^'^^<^) (164) 

and thus we obtain for the wave numbers /c„ 

kn = -^n, nGZ\{0}, (165) 
In (a) 

wherein we have taken into account that A = A;q = is not an eigenvalue for the 
Dirichlet case. In contrast to the Dirichlet case, A = /cq = is an eigenvalue for 
Neumann boundary conditions at both interval ends, and the nonzero wave numbers 
coincide with the Dirichlet case (11650 . For Robin boundary conditions at both interval 
ends, we get 

F{k) = l- ^P'-'^}}P'-% ^^^<i). (166) 
^ ' (pi + ik) ip2 + ik) 
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Again, by a comparison of f ll65p with (I162p . we conclude for the Maslov indices for pure 
Dirichlet and Neumann boundary conditions 

fin = for neN (167) 

and additionally /iq = for the Neumann case. For the Robin boundary conditions on 
both interval ends the Maslov indices have to be individually calculated for each n G Nq 
by (fT66D . 

The above considerations underline the fact that the form of the S-matrix fll55p 
corresponds to a pure ring system as in the case of the negative Laplacian —A, see 
|44] . The occurrence of possible noninteger "Maslov indices" originates simply from the 
fact that we have a discontinuous crossover by turning once around in the ring system 
(one- dimensional torus) in contrast to the "usual" continuity requirement of the wave 
function, see e.g. [57] . 



Example 16.3 We shall present a trace formula for the time- evolution operator U{t) 
in l[21\) and [2^] for H^k acting on a single edge with the assigned interval I = [1, b]: 

oo 

(t/(t)0) (x) := M^) {^n, 0) e-''^"*, G L\I, dx), (168) 

n=— oo 

with the eigenvalues kn = ^{n-\-c), n G Z, c G [0, 1), and the normalized eigenfunctions 

^„(x) = ^^e^'^"'"^ with neZ. (169) 
\/x\nb 

For the corresponding (not retarded) integral kernel of U{t) we get by \2^^\ p. 20] 
(in a distributional sense acting on T>{I) C L'^{I,dx) identified by the continuous 
representatives ; g{x, Xq] t) := [Inx — Inxo — t]) 

°° _ „icg(x,xo;t) °° 

K{x,xo;t):= J2 Mx)i^nMe-''"' = e^^^^'^"^*)" 

n=— oo ^ n=—oo 

Y 5ig{x,xo;t) + 2Tm) (170) 



In b ^/xxo 

c 

n=— oo 

If we take the trace ofU{t), we obtain with U70\) (by defining the "period" T = ln6 [see 
^] and the Maslov index fi := Ac) 

b 

oo 



n=— oo 

oo 



K{x,x;t)dx = Y ^~ 
1 "=-°° (171) 

oo 
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// we now choose a test function h of Hr ( definition \1 5. 1\) with an arbitrary r > 0, we 
get by the identity U71\) and the symmetry of h the trace formula 

oo 

/oo 
h{t)TTU{t)dt= J2 ^(^^) 
n=— oo 

— oo 

oo 

= Th{0) +tJ2 (e"'5/^"/i(nT) + e'5^"/i(-nT)) (1^2) 



n=l 



cos (^■^A*'^) h{nT). 

n=l 

We recall that the S -matrix for this quantum graph is S{A,B) = e'"^^^^ = e~^^^ (see 
Ij^TJ^), and the length of the (single) edge is I = £, = Irij = Info. Since we have 
a directed edge, there is only one possibility for the orientation of the periodic orbits 
and, therefore, the periodic orbits can be labelled by the natural numbers, and the 
corresponding lengths of the periodic orbits are [„ = rain 6 and all are multiples of one 
primitive periodic orbit with length li = ln6 = T. For the amplitude functions we get 
(see [37\) An = lie^^™ = Te-^f'^". Applying DM) we get HW . which confirms the 
trace formula in theorem \15.2[ 

Example 16.4 We shall present a trace formula for the kernel K{x, xq', t) of the unitary 
evolution operator e^'*^BK of -^bk Dirichlet boundary conditions (D) on a single 
edge e with assigned interval I = [l,b]. The eigenvalues are given by U65\) . thus the 
(Feynman-)kernel reads 

oo 

K{x, xo; t) := Y,M^)i^ni^o)e-"'"' (173) 

n=l 

with the normalized eigenf unctions 

^„(x):=W- ^eN, (174) 



where \ := Info denotes the length of the edge e. Using a suitable addition theorem 
for trigonometric functions, we get two alternative expressions for K{x,Xo',t) (see 
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K{x,xo;t) 



1 



■7p 



[2 

7p 



7p 




gi2n7rg^p I 1 



. (inx — Inxo + ni^p)^ 



At 



^gi(2n+i)^ exp I 1 



(inx + Inxo + n\^^)' 



:i75) 



+ 6^2'^" exp I 1 



. (in X — In xq — nl^p) ' 



At 



^gi(2n+i)^ exp I 1 



. (Inx + Inxo - n\^y 



At 



where Q^{z,t) denotes the Jacobi theta function and we have defined 

1 
2 

1 for n > 0. 



for n = 



(176) 



:= 2[ = 2 In 6 is the length of the primitive periodic orbit 7^ of the corresponding 
classical system. Notice that the summands in the second identity in U75\) can be 
interpreted as contributions of free particle kernels at a fixed time t corresponding to 
the four types of paths p{xo,x) (see sectionl^ joining xq and x (see e.g. [5M l5E\). For 
this reason, we define (cf. and [23 p. 30]) 



Kp(x, Xq] t) :- 



2Vi7rt 



exp I 1 



At 



(177) 



where ip{xQ,x) is the length of the path p{xq,x) (see ^^), and we then get by Iil75\) 
K{x,xo;t) = ^ exp {innp^^^^^)) Kp{x,xo;t) 



/XXq 

1 

/XXo 



p{xo,x) 

E 

p{xo,x) 



exp 



• '^f^p(xo,x) 



Kp(x, Xq] t), 



(178) 



where the sum comprises all possible paths p{xo,x) joining xq and x, and np(xo,x) is 
defined as the number of reflections of the path p{xo,x) at the "hard wall" interval 
endpoints 1 and Info. Hp(xo,x) denotes the Maslov index of the path p{xo,x) which is 
given by fip(xo,x) = 2np(2;o^a;) mod4 in agreement with the "usual" Maslov index for the 
one- dimensional billiard system corresponding to the negative Laplacian (see [5^j, and 
with M61\) (in lil61/\] the Maslov index corresponds to periodic orbits). 

Example 16.5 Finally, we shall present an explicit trace formula (heat kernel) for a 
single edge with assigned interval I = [1,6] for with Dirichlet boundary conditions 
{D). We calculate the trace of the heat kernel of e~^^^^ (replacing t by —it in U73\) ) 



Khix, xo] t) := K{x, xo] -it) = ^ tJj{x)^^{xo)e 



■ k„ t 



(179) 
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directly and then compare the result with the trace formula U37\) . Therefore, we recall 
that the wave numbers of H^^^ with {D) are explicitly given by U65\) . Thus, we obtain 
for the trace of the heat kernel (setting £, := := Info and \ 54,P- 371] j 

b 



(180) 




Notice that the sums in U80\) are absolutely convergent whereas in U71\) the sums are 
convergent in the topology o/P'(M) (in a distributional sense). In order to compare this 
result with USI]) . we recall that C{1, In 6; k') and the S-matrix S'p) for the Dirichlet case 
is given by (see U21\) and U63\) ) 

Infc ^jr 

^(% = -l2x2 and CilMb;k')= \ 1. (181) 



2^+lnb 



It is a simple calculation that the multiplicity qq of the eigenvalue one of S"jj^C{l, In b; k') 
is go = for any A;' G R \ {0}. Furthermore, it is obvious that the order N of the zero 
with wave number kQ = of F{k) in M18) is N = 1, thus Qq — = — |. The 
multiplicities of the wave numbers k^ are Qn = I for n G N. Since a Dirichlet boundary 
condition corresponds to the classical "hard wall" boundary condition, we conclude that 
the periodic orbits 7 are given by all multiples of one primitive periodic orbit 7^ with 
primitive periodic orbit length l^^ = 2 In 6. For the amplitude functions in M31\] 
we obtain = (see \3l\). Furthermore, it is obvious that ImSo = 0. Using the 
test function h{k) := * we obtain the Fourier transform h{x) = ^^^^^47. Inserting 
these quantities in the trace formula I^137\) we get the trace formula I^180\) . which again 
confirms the trace formula UST]) . We remark that from the small-t asymptotics U80\) 
one obtains directly the Weyl asymptotics I^138i) using a proper Karamata-Tauberian 
theorem. 



17. Summary and conclusions 

We have studied tlie quantization of tlie extraordinarily simple classical Hamiltonian 
Helix, p) = xp about which Berry and Keating [I3l[l4] speculated that some quantization 
of it might yield the hypothetical Hilbert-Polya operator [11 |2l [3l HJ [TOl [TTl [12] possessing 
as eigenvalues the nontrivial Riemann zeros. Two quantum Hamiltonians respectively 
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Schrodinger operators have been considered: The original Berry-Keating operator 
i^BK '■= (^3^ + I) obtained from Hd by Weyl ordering, and the so-called "squared" 
Berry-Keating operator -f^BK '■— ~ ^-^Hx ~ i which is a special case of the famous 

Black-Scholes operator [33| IM] used in the financial theory of option pricing. 

In section [21 we have given a complete description of the quantum dynamics 
generated by i^BK acting in the Hilbert space L^(M>,da;). While the one-dimensional 
quantum system governed by //bk possesses many interesting properties, one of our main 
results of section |2] is that the spectrum of H^k is purely continuous corresponding to 
scattering states. Since there are no bound states corresponding to a discrete spectrum, 
it is obvious that this specific quantization of the Berry-Keating operator cannot possess 
the Riemann zeros as part of its spectrum. Let us point out, however, that in the simple 
example 116. II we have studied the quantum dynamics of Hbk for the particular square- 
integrable wave function ( I14ip for which it turns out that the spectral decomposition 
consists of a continuous wave number spectrum, A; G M, into which there are embedded 
infinitely many absorption lines located exactly at the wave numbers /c„ = r„ G M 
corresponding to the nontrivial Riemann zeros satisfying the Riemann hypothesis. To 
our knowledge, there is, however, no relation of this occurrence of the Riemann zeros to 
the absorption spectrum interpretation of Connes [TT| [T2]. 

Analogous results have been obtained in section H] for the "squared" Berry-Keating 
operator -f^BK acting in the Hilbert space L^(M>, dx). We have proved that in this case 
the spectrum is purely continuous too and thus there holds again a "no-go theorem" 
with respect to the identification of with the hypothetical Hilbert-Polya operator. 

In the main part of our paper, we have dealt with the quantum dynamics of -f^BK 
respectively -^bk compact quantum graphs introduced in section [6l After having 
defined the Berry-Keating operator i^BK on compact graphs in section [TJ we have 
given in proposition 18.11 a complete classification of all self-adjoint extensions of i^BK on 
compact quantum graphs in terms of two matrices A and B satisfying the conditions 
(168|) . In proposition 19 . 1 1 we have established a secular equation valid for any self-adjoint 
realization in form of a determinant whose zeros determine the discrete spectrum of 

In the sections [T0l - IT3l we have studied the quantization of the "squared" Berry- 
Keating operator i^BK compact quantum graphs. Proposition 111.11 provides the 
complete classification of all self-adjoint realizations of -f^BK again in terms of matrices 
A and B satisfying in this case the conditions (1971) . For the discrete spectrum of i^BK^ 
we have given in proposition 112.11 the corresponding secular equation for X = k"^ 0. 
The zero eigenvalue A = A;^ = of i^BK plays a special role which we have discussed in 
section [13] leading to the additional secular equation fll2ip . Furthermore, in section [T^ 
we have discussed the conditions under which -f^BK square of i^BK in the sense of 

(IHTl) . see proposition 114.11 

Based on the results derived in the previous sections, we have been able in section 
[TSlto state several theorems. In the theorems 115. 21 and 115. 31 we have given an exact trace 
formula for H^k respectively i^BK ^ large class of test functions h{k) belonging to 
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the space Hj. defined in definition 115. II The trace formulae establish a deep connection 
between the eigenvalue spectra of i^BK respectively -f^BK ^"^^ length spectra of the 
periodic orbits of the corresponding classical dynamics. 

As an important consequence of the trace formulae, we have derived the Weyl 
asymptotics for Hbk (for -^^bk have used results in [37J). The Weyl asymptotics 
of these operators have been given in the theorems 115.51 and 115. 4[ respectively. A 
comparison with the expected Weyl asymptotics (HHj) respectively (l32l) for the nontrivial 
Riemann zeros demonstrates clearly that neither Hbk nor i^BK "^^^ yield as eigenvalues 
the nontrivial zeros of the Riemann zeta function if these operators are self-adjoint 
realizations on any compact quantum graph, see theorem 115.61 

Finally, we have presented in section [16] four simple examples illustrating some 
aspects of the quantum dynamics of Hbk respectively H^^. 



Acknowledgements 

We would like to thank Wolfgang Arendt, Jens Bolte and Jan Eric Strang for very 
helpful discussions and fruitful hints and Wolfgang Arendt for drawing our attention 
to [271 [23 [29]. S. E. would like to thank the graduate school "Analysis of evolution, 
information and complexity" of the Land Baden- Wiirttemberg for the stipend which has 
enabled this paper. 

[1] In the comment (included in his collected papers [1]) to his paper on the explicit formulae of 
number theory [3], A. Weil writes that E. Hellinger once told him the following: At the occasion 
of his first course on his theory of integral equations, Hilbert said after having demonstrated 
that the characteristic values (eigenvalues) of a symmetric kernel are real: "Et avec ce theoreme, 
Messieurs, nous demontrerons I'hypothese de Riemann". Since Hellinger was a student of Hilbert 
in Gottingen (since 1904) and later until 1909 Hilbert's assistant, we can assume that the cited 
statement by Hilbert dates back to these years. 

[2] A. Weil. Commentaire 'Sur les "formules explicites" de la theorie des nombres premiers'. In 
Scientific works. Collected papers., volume II (1951-1964), page 527. Springer, New York, 1979. 

[3] A. Weil. Sur les "formules explicites" de la theorie des nombres premiers. Comm. Sem. Math. 
Univ. Lund [Medd. Lunds Univ. Mat. Sem.], 1952:252-265, 1952. 

[4] See Odlyzko's correspondence with Polya: |http : //www . dtc . umn . edu/~odlyzko/polya/ index .html! 

[5] A. Selberg. Harmonic analysis and discontinuous groups in weakly symmetric Riemannian spaces 
with applications to Dirichlet series. J. Indian Math,. Soc. (N.S.), 20:47-87, 1956. 

[6] J. Hadamard. Les surfaces a courbures opposees et leurs lignes geodesiques. Journal de 
Mathematique pure et appliquees, 4:27-73, 1898. 

[7] J. Hadamard. Sur le billard non-Euclidien. Soc. Sci. Bordeaux, Proc. Verbaux, 147, 1898. 

[8] M. C. Gutzwiller. Classical quantization of a Hamiltonian with ergodic behaviour. Phys. Rev. 
Lett, 45:150-153, 1980. 

[9] M. C. Gutzwiller. Periodic orbits and classical quantization conditions. J. Math. Phys., 12:343- 
358, 1971. 

[10] M. V. Berry. Riemann's zeta function: a model for quantum chaos? In T. H. Seligman and 
H. Nishioka, editors. Quantum Chaos and Statistical Nuclear Physics, volume 263 of Lecture 
Notes in Physics, pages 1-17. Springer, New York, 1986. 

[11] A. Connes. Formule de trace en geometric non-commutative et hypothese de Riemann. C. R. 
Acad. Sci. Paris Ser. I Math.. 323:1231-1236, 1996. 



The Berry- Keating operator on L (M>, dx) and on compact graphs 38 
[12 
[13 



[14 

[is: 

[16 
[17] 

[18; 

[19 
[20 
[21 

[22 

[23 

[24 

[25 
[26 

[27 

[28 

[29 

[30 

[31 

[32 

[33 

[34 
[35 

[36 

[37] 



A. Connes. Trace formula in noncommutative geometry and the zeros of the Riemann zeta function. 

Selecta Math. (N.S.), 5(1):29-106, 1999. 
M. V. Berry and J. P. Keating. H = xp and the Riemann zeros. In I. V. Lerner, J. P. Keating, and 
D. E. Khmelnitskii, editors, Supersymmetry and Trace Formulae: Chaos and Disorder, volume 
263, pages 355-367. Plenum Press, New York, 1999. 
M. V. Berry and J. P. Keating. The Riemann zeros and eigenvalue asymptotics. SIAM Rev., 
41(2):236-266, 1999. 

T. Kottos and U. Smilansky. Quantum chaos on graphs. Phys. Rev. Lett., 79(24):4794-4797. 
T. Kottos and U. Smilansky. Periodic orbit theory and spectral statistics for quantum graphs. 

Ann. Phys. (NY), 274:76-124, 1999. 
J.-P. Roth. Spectre du laplacien sur un graphe. C. R. Acad. Sci. Paris Ser. I Math., 296(19):793- 

795, 1983. 

J. von Below. A characteristic equation associated to an eigenvalue problem on -networks. 

Linear Algebra AppL, 71:309-325, 1985. 
P. Exner and P. Seba. Free quantum motion on a branching graph. Reports on Mathematical 

Physics, 28:7-26, 1989. 

P. Exner and O. Post. Convergence of spectra of graph- like thin manifolds. J. Geom. Phys., 
54(1):77-115, 2005. 

0. Post. Spectral convergence of quasi-one-dimensional spaces. Ann. Henri Poincare, 7(5):933- 
973, 2006. 

R. Carlson. Linear network models related to blood flow. In Quantum graphs and their 
applications, volume 415 of Contemp. Math., pages 65-80. Amer. Math. Soc, Providence, RI, 
2006. 

1. M. Gelfand and G. E. Schilow. Verallgemeinerte Funktionen (Distributionen). I: 
Verallgemeinerte Funktionen und das Rechnen mit ihnen. Hochschulbiicher fiir Mathematik, 
Bd. 47. VEB Deutscher Verlag der Wissenschaften, Berlin, 1960. 

M. Reed and B. Simon. Methods of Modern Mathematical Physics. IL Fourier Analysis, Self- 
Adjointness. Academic Press, New York, 1975. 

E. C. Titchmarsh. Introduction to the Theory of Fourier Integrals. Oxford, 1948. 

J. Weidmann. Lineare Operatoren in Hilbertrdumen. Teil II. Mathematische Leitfaden. 

[Mathematical Textbooks]. B. G. Teubner, Stuttgart, 2003. 
W. Arendt. Gaussian estimates and interpolation of the spectrum in L^. Differential Integral 

Equations, 7(5-6):1153-1168, 1994. 
W. Arendt. Spectrum and growth of positive semigroups. In Evolution equaMons (Baton Rouge, 

LA, 1992), volume 168 of Lecture Notes in Pure and Appl. Math., pages 21-28. Dekker, New 

York, 1995. 

W. Arendt and B. de Pagter. Spectrum and asymptotics of the Black-Scholes partial differential 
equation in (L\ L°°)-interpolation spaces. Pacific J. Math., 202(l):l-36, 2002. 

C. Grosche and F. Steiner. Handbook of Feynman path integrals, volume 145 of Springer Tracts 
in Modern Physics. Springer- Verlag, Berlin, 1998. 

S. Endres and F. Steiner. A simple infinite quantum graph. Ulmer Seminare 2009, 14:187-201, 
2009. 

M. C. Gutzwiller. Chaos in classical and quantum mechanics, volume 1 of Interdisciplinary Applied 
Mathematics. Springer- Verlag, New York, 1990. 

F. Black and M. S. Scholes. The pricing of options and corporate liabilities. Journal of Political 
Economy, 81(3):637-54, 1973. 

R. C: Merton. Theory of rational option pricing. Bell Journal of Economics, 4(1):141-183, 1973. 
A. Terras. Harmonic analysis on symmetric spaces and applications. I. Springer- Verlag, New 
York, 1985. 

V. Kostrykin and R. Schrader. Kirchhoff's rule for quantum wires. J. Phys. A, 32:595-630, 1999. 
J. Bolte and S. Endres. The trace formula for quantum graphs with general self adjoint boundary 



e Berry-Keating operator on (M>, dx) and on compact graphs 



39 



conditions. Annales Poincare, 10(1): 189-223, 2009. 
[38] J. Weidmann. Lineare Operatoren in Hilbertrdumen. Teil 1. Mathematische Leitfaden. B. G. 
Teubner, Stuttgart, 2000. 

[39] M. Reed and B. Simon. Methods of modem mathematical physics. I. Academic Press Inc., New 

York, second edition, 1980. Functional analysis. 
[40] P. Kuchment. Quantum graphs. I. Some basic structures. Waves Random Media, 14:S107-S128, 
2004. 

[41] M. Reed and B. Simon. Methods of Modem Mathematical Physics. IV. Analysis of Operators. 

Academic Press, New York, 1978. 
[42] W. N. Everitt and L. Markus. Complex symplectic spaces and boundary value problems. Bull. 

Amer. Math. Soc. (N.S.), 42(4):461-500, 2005. 
[43] V. Kostrykin and R. Schrader. Laplacians on metric graphs: eigenvalues, resolvents and 

semigroups. In G. Berkolaiko, R. Carlson, S. A. Fulling, and P. Kuchment, editors, Quantum 

graphs and their applications, volume 415 of Contemp. Math., pages 201-225. Amer. Math. Soc, 

Providence, RI, 2006. 

[44] V. Kostrykin and R. Schrader. The inverse scattering problem for metric graphs and the traveling 

salesman problem, preprint, arxiv:math-ph/0603010, 2006. 
[45] R. Carlson. Inverse eigenvalue problems on directed graphs. Trans. Amer. Math. Soc, 

351(10):4069-4088, 1999. 

[46] B. Thaller. The Dirac equation. Texts and Monographs in Physics. Springer- Verlag, Berlin, 1992. 
[47] S. A. Fulling, P. Kuchment, and J. H. Wilson. Index theorems for quantum graphs. J. Phys. A, 
40:14165-14180, 2007. 

[48] E. Schrodinger. A method of determining quantum-mechanical eigenvalues and eigenfunctions. 

Proc. Roy. Irish Acad. Sect. A., 46:9-16, 1940. 
[49] L. Infeld and T. E. Hull. The factorization method. Rev. Modern Physics, 23:21-68, 1951. 
[50] P. Kurasov and M. Nowaczyk. Inverse spectral problem for quantum graphs. J. Phys. A, 38:4901- 

4915, 2005. 

[51] V. Kostrykin and R. Schrader. Quantum wires with magnetic fluxes. Comm. Math. Phys., 
237:161-179, 2003. 

[52] V. Kostrykin, ,J. PotthofF, and R. Schrader. Heat kernels on metric graphs and a trace formula. 
In F. Germinet and P.D. Hislop, editors, Adventure in Mathematical Physics, volume 447 of 
Contemp. Math., pages 175-198. Amer. Math. Soc, Providence, RI, 2007. 

[53] J. Karamata. Neuer Beweis und Verallgemeinerung einiger Tauberian-Satze. Math. Z., 33:294- 

299, 1931. 

[54] W. Magnus, F. Oberhettinger, and R. P. Soni. Formulas and theorems for the special functions 
of mathematical physics. Third enlarged edition. Die Grundlehren der mathematischen 
Wissenschaften, Band 52. Springer, New York, 1966. 

[55] A. Einstein. Zum Qantensatz von Sommerfeld und Epstein (in German). Verh. Deutsche Phys. 
Ges., 19:82-92, 1917. 

[56] J. B. Keller. Corrected Bohr-Sommerfeld quantum conditions for nonseparable systems. Ann. 

Physics, 4:180-188, 1958. 

[57] H.-J. Stockmann. Quantum chaos. Cambridge University Press, Cambridge, 1999. An 
introduction. 

[58] F. Barra and P. Gaspard. Transport and dynamics on open quantum graphs. Phys. Rev. E (3), 
65(1, part 2):016205, 21, 2002. 



